Coherent population transfer in a chain of tunnel coupled quantum dots* 
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We consider the dynamics of a single electron in a chain of tunnel coupled quantum dots, exploring 
the formal analogies of this system with some of the laser-driven multilevel atomic or molecular 
systems studied by Bruce W. Shore and collaborators over the last 30 years. In particular, we 
describe two regimes for achieving complete coherent transfer of population in such a multistate 
system. In the first regime, by carefully arranging the coupling strengths, the flow of population 
between the states of the system can be made periodic in time. In the second regime, by employing a 
"counterintuitive" sequence of couplings, the coherent population trapping eigenstate of the system 
can be rotated from the initial to the final desired state, which is an equivalent of the STIRAP 
technique for atoms or molecules. Our results may be useful in future quantum computation schemes. 

PACS numbers: 03.67.-a, 73.63.Kv, 73.23.Hk 



I. INTRODUCTION 

Population transfer in multistate quantum systems has 
been an active topic of research over the last half a cen- 
tury. In the context of atomic and molecular physics, co- 
herent population transfer in optically-driven multilevel 
systems has been studied since the invention of lasers 
1]. Usually, the objective is to transfer the population 
from the initial to a well defined final state of the atom 
or molecule, via one or more intermediate states, while 
minimizing the loss of population through or its accu- 
mulation on the intermediate states. In early theoretical 
work, Shore and collaborators have studied population 
transfer in multilevel systems driven by resonant laser 
fields @, 0j IH- In particular, they have found that it is 
possible to arrange the coupling strengths between the 
adjacent states in such a way that the system becomes 
analogous to a spin- J in a magnetic field, whose dynamic 
evolution is known to be periodic for any J [3| . This cou- 
pling scheme was therefore named spin-coupling. 

Later, Hioe, Eberly, Bergmann and collaborators dis- 
covered the technique of stimulated Raman adiabatic 
passage (STIRAP) for three- level atomic/molecular sys- 
tems Q. They have identified a specific eigenstate of 
the system, the so-called coherent population trapping 
(CPT) state, which contains a superposition of the ini- 
tial and final states, and dates back to Alzetta et al. 
and Arimondo and Orriols The STIRAP technique 
is then based on first preparing the system in its initial 
bare state, which coincides with the CPT state, and then 
adiabatically rotating the CPT state towards the desired 
final bare state of the system. This techniques has been 
subsequently polished Q and extended to multilevel sys- 
tems 0, 0, Il0j | with the active participation of Bruce W. 
Shore. 

While the above studies were conducted in the con- 
text of multilevel atoms or molecules, here we show 
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that similar effects can be found in the context of quan- 
tum transport in arrays of tunnel-coupled quantum dots 
[Til . H3 , Ha , Often referred to as artificial atoms, 

semiconductor quantum dots offer an unprecedented pos- 
sibility of constructing at will and exploring situations 
ranging from practically single atom to a fully solid state 
many-body systems [15| . The nanofabrication possibili- 
ties of tailoring structures to desired geometries and spec- 
ifications, and controlling the number and mobility of 
electrons confined within a region of space, are some of 
the features that make these structures unique tools for 
the study of a variety of preselected set of phenomena, 
including the coherent population transfer in multistate 
systems. 

Given the controllable quantum properties of the elec- 
trons in such structures, the possibility of their applica- 
tion to schemes of quantum computers (QCs) [16| has 
not escaped attention [Tt], EH E|- The qubits of the 
QD-array based QC would be represented by the spin- 
states of single electrons confined in individual QDs, 
with the two-qubit nearest-neighbor coupling mediated 
by the controlled spin-exchange interaction [l7], [H] . One 
of the main difficulties with the existing proposals for 
integrated solid-state based QCs is that there is no effi- 
cient way of transferring the information between distant 
qubits. We consider here a single-electron tunneling in 
a one-dimensional array of QDs and establish the con- 
ditions under which the complete transfer of the elec- 
tron wavepacket between two distant locations can be 
achieved. Our findings could therefore be relevant to the 
reliable information exchange between distant parts of an 
integrated quantum computer [20| . 

In Section [II] we outline the mathematical formalism 
describing a chain of QDs, in terms of which, in Section 
IIII1 we present the theory of coherent propagation and 
periodic oscillations of the electron wavepacket between 
the two ends of the chain. The single-electron transfer 
via an equivalent of multistate STIRAP is discussed in 
Section HVl In Section [V] we describe an envisioned im- 
plementation of a scalable quantum computer, followed 
by the concluding remarks. 
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FIG. 1: Schematic drawing of the chain of tunnel-coupled 
QDs. 



II. MATHEMATICAL FORMALISM 

We consider electron transport in a linear array of 
N nearly identical QDs which are electrostatically de- 
fined in a two-dimensional electron gas by means of 
metallic gates on top of a semiconductor heterostructure 
(GaAs/AlGaAs) [H Gjjj]. This system is described by 
the extended Mott-Hubbard Hamiltonian [12, EH , which 
in its most general form is given by 

H = Y2 £ joca] a a ja + - J^Unjinj - 1) 



3,<x 



i<3 



where at and 

Qjcx axe the creation and annihilation 
operators for an electron in state a with the single- 
particle energy Ej a , U is the on-site Coulomb repulsion, 
rij = dj a aj a the total electron number operator of 
the jth dot, tij^ a are the coherent tunnel matrix elements 
between dots i and j, and Vij is the interdot electrostatic 
interaction. In general, the index a refers to both or- 
bital and spin states of an electron. In the tight-binding 
regime, when the on-site Coulomb repulsion and single- 
particle level-spacing Ae are much larger than the tun- 
neling rates, U > Ae 3> Uj ja , only the equivalent states 
of the neighboring dots are tunnel- coupled to each other 
pll |. In the absence of a magnetic field, we can thus limit 
our consideration only to a single doubly- (spin-) degen- 
erate level per dot (a G {Ti.I})j assuming further that 
the tunneling rates do not depend on the electron spin. 

In this paper we are concerned with single-electron dy- 
namics, considering a situation in which a preselected QD 
is initially doped with one mobile electron, while all of the 
other dots of the chain are empty, as indicated in Fig. [T] 
Our aim is to determine the conditions under which the 
complete coherent transfer of the electron between the 
two ends of the chain can be achieved. The population 
transfer in this system is mediated by the tunneling be- 
tween the neighboring QDs. The individual tunneling 
rates tj = fjj+i are determined by the voltages applied 
to the gates defining the corresponding interdot tunneling 
barriers. A chain of N tunnel-coupled QDs doped with a 
single electron is described by the following Hamiltonian, 

H lc = ^e J al a a ja +J2 t j( a l a a 3+i,c l + a ja a' < j+la ), (2) 



which obviously does not contain terms responsible for 
electrostatic interactions. Since this Hamiltonian pre- 



serves the electron number and its spin, the total state- 
vector of the system reads 



Wr)>=5>jM \j a ), 



(3) 



where \j a ] 



= „t 



|0i,...,0jv) denotes the state with one 



electron having spin a at the jth dot. The time-evolution 
of the system is governed by the Schrodinger equation 
i \ip) — H\ c \ip) (S- = 1), which yields 



e J A j 



(4) 



where to = — 0. Obviously, the two sets of these am- 
plitude equations with a =] and a ={ are equivalent and 
decoupled from each other. As a result, if the electron 
is prepared in an arbitrary superposition of spin up and 
spin down states, \ip) = A| + A 1 - the two parts 
of the wavefunction evolve symmetrically and indepen- 
dently of each other. This assertion is valid as long as 
all the uncontrollable spin-flip processes are vanishingly 
small on the time scale of t^ 1 . In semiconductor QDs, the 
spin decoherence originates mainly from the spin-phonon 
coupling, as well as the coupling of the electron spin with 
the nuclear spins of the surrounding crystal (hyperfine 
interaction) or stray magnetic fields. The first decoher- 
ence mechanism is suppressed at low temperatures 211. 
at which the density of crystal phonons is negligible 22 1 . 
As for the uncontrollable hyperfine interactions, exper- 
imental measurements indicate spin-relaxation times in 
excess of 100 //s, which can be further improved by ap- 
plying moderate magnetic fields or polarizing the nuclear 
spins [23[. Another mechanism for decoherence in the 
process of electron (charge) transfer in our system origi- 
nates from the structure imperfections and gate voltage 
fluctuations, which cause uncertainty in the intradot en- 
ergy levels and interdot couplings. These fluctuations, 
however, are typically slow on the time scale of f _1 , and 
the resulting disorder in the system may be considered 
frozen during its dynamic evolution, as we have discussed 
in a previous publication [20l |. 

Let us write the Hamiltonian for the electron with spin 
a in the matrix form 



H? e = 



ei h o 

tl £2 t2 

h e 3 



£JV- 
tN- 



en 



(5) 



which is obviously tridiagonal. Inspection of the ampli- 
tude equations 0$ or the Hamiltonian ([5]) indeed verifies 
that our system is formally analogous to the laser-driven 
multilevel atomic or molecular systems studied by Shore 
and coworkers @, 0, [j] and Bergmann, Shore and others 
[EH II, SOU- Here, the tunneling rates tj between states 
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| j) and \j + 1) play the same role as the Rabi frequen- 
cies of the laser fields acting on the atomic transitions 
li) \.i + l)i while the energies Sj of states |j) corre- 
spond to the cumulative detunings of the atomic levels. 
In the following Sections, we describe two methods for 
achieving complete population transfer from the initial 
|1) to the final \N) state of the system, which turn out 
to be the counterpart of those in Refs. [|[ and [Toj . 



III. PERIODIC OSCILLATIONS OF 
POPULATION BETWEEN THE TWO END 
STATES 

In this Section we consider the electron wavepacket 
dynamics in the chain with static couplings between the 
dots. Assume that at time r = the electron is local- 
ized on the first dot, ^"(O)) = |l a ), and the tunnel 
couplings are switched on. This switching should be fast 
enough on the time scale of t _1 , so that no appreciable 
change in the initial state of the system occurs during 
the switching time r sw , but slow on the time scale of 
so that no nonresonant coupling between the dots is in- 
duced: e -1 < t sw < t . The aim is to determine the set 
of couplings between the states of the systems which will 
achieve a complete transfer of the electron population 
from the initial to the final dot. 

To determine the time-evolution of the state vector 
(|3|) we need to solve the eigenvalue problem Hf e \ip a ) — 
\\4> a ) which will yield the eigenvalues Afc and corre- 
sponding eigenvectors of the Hamiltonian ([5]). The 
state vector \4> a {T)) at any time r > is given by 



N 



N 



k j 

(6) 

Note that the matrix in Eq. ([5]) has the form of the 
tridiagonal Jacobi matrix. It is natural to first con- 
sider the case of equal tunneling rates between the dots: 
tj = t. Assuming equal energies Sj = e and making 
the transformation A" — > A^e %eT ', which is equivalent 
to the interaction picture, we find that the determinant 
I? at (A) = det(iJf c — AI) is identical to the Chebyshev 
polynomial of the second kind, which can be expressed as 
f jv(A) = n^ =1 (A — Afc). The eigenenergies of the system 
are then given by the roots of this polynomial, namely 



An 



2t cos 



kn 
N + l 



while the corresponding eigenvectors are 

N 



N 



TiE sin G^ri) ^ 



Using Eq. © and the initial conditions A\ — 1 and 
Aj = for j = 2, 3, . . . N, we obtain the solutions for the 



amplitudes as, 



N- 



x sin 



N 

iE ex p 

fe=i 

N + l 



-i2tr cos 



kn 



kit 
N+l 



N+l 



(7) 



It is thus evident that the eigenstates of the coupled 
system oscillate with incommensurate frequencies corre- 
sponding to the roots Afc of I? at, which in fact become in- 
creasingly densely spaced with increasing N. As a conse- 
quence, the system never revives fully to its initial state, 
as is illustrated in Fig.[2Ua). 

Clearly, it is highly desirable to tailor the parameters 
of the system so as to achieve a non-dispersive transfer 
of the single-electron wavepacket between the two ends 
of the chain. Recall from the theory of angular momen- 
tum that a spin-J particle subject to a constant mag- 
netic field exhibits Larmor precession about the field di- 
rection. In particular, if one chooses the quantization 
direction along an axis perpendicular to the magnetic 
field direction and prepares the particle in its lowest spin 
eigenstate \J,M = — J), it will oscillate between this 
initial and the final state \J,M = J) in a perfectly 
periodic way. The matrix elements for the transitions 
| J, M) «-> |J, M + 1) between the neighboring states are 
proportional to y/{J — M)( J + M + 1). It is therefore 
clear that with the appropriate choice of the interdot 
tunneling matrix elements, the dynamics of the single- 
electron in a chain of QDs can mimic that of a spin-J in 
a magnetic field. Indeed, if we formally set N = 2 J + 1 
and j = J + M + 1, the tunneling rates tj should be ar- 
ranged according to tj = ty/ (N — j)j for j = 1, N— 1. 
Then again, by exploring the properties of the Jacobi 
polynomials, we find equally spaced eigenenergies of the 
system, 

Afc =t(2k- N- 1), 

while the corresponding eigenvectors can be expressed 
through the rotation matrices commonly used in the rep- 
resentation theory of angular momentum. With the ini- 
tial conditions A\ — 1 and Aj — for j — 2,3, ... N, for 
the amplitudes of the state- vector ([3]), we then obtain 
simple analytic expressions given by the binomial form 



-4'? 



N - 1 
3-1 



1/2 



-isin (tr)^'^ cos (tr) 



(8) 



Since the eigenstates of the system have commensurate 
energies Afc, the electron wavepacket oscillates in a per- 
fectly periodic way between the first and the last dots, 
whose occupation probabilities are given, respectively, 
by \A?\ 2 = cos(tr) 2(JV-1) and l^l 2 = sin (*r) 2(JV-1) , 
which is illustrated in Fig.[2jb). In particular, if at time 
t = 7i"/ (2t) the tunneling rates are suddenly switched off, 
we obtain \Af\ 2 = and lA^I 2 = 1, i.e. complete pop- 
ulation transfer from the initial to the final state of the 



FIG. 2: Time-evolution of a single-electron wavepacket in a chain of N = 9 QDs with static tunneling rates, (a) Population flow 
in the chain with equal interdot tunneling rates tj = t (shown in the inset), (b) Population flow in the chain with spin-model 
tunneling rates tj = t\J (N — j)j (shown in the inset). The time r is in units of t . 



system. In a somewhat abstract sense, the behavior of 
the system is thus similar to that of a two-level system 
subject to a tt pulse. Let us note at this point that the 
population transfer between the two ends of the chain 
can be achieved most straightforwardly by sequentially 
pulsing the tunneling rates between the first and second 
dots for time tj = ir/(2t 1 ), then the second and third 
dots for time t 2 = tt/ (2t 2 ), etc till reaching the Nth. dot, 
which is equivalent to applying a sequence of tt pulses 
in a multistate atomic system. In the scheme described 
above, however, all the interdot tunnelings are switched 
on and then off simultaneously, realizing thereby a fast 
and efficient transfer of the electron from the first to the 
last QD. 



IV. ADIABATIC POPULATION TRANSFER 
BETWEEN THE TWO END STATES 

While the above tunneling schemes, involving a se- 
quence of 7T pulses or an effective collective tt pulse, re- 
quire both, careful control of the individual tunneling 
rates and their timing, in this Section we describe a ro- 
bust adiabatic method for population transfer which is 
not very sensitive to small uncertainties in the interdot 
tunneling rates. Recall that a three- level atom inter- 
acting with two laser fields, under the condition of two- 
photon (Raman) resonance, possesses a coherent popula- 
tion trapping (CPT) state, which is decoupled from both 
laser fields [7[. Equivalently, for a chain of three tunnel- 
coupled quantum dots, assuming equal energies Sj = e, 
the eigenstate of Hamiltonian ([5]) with zero eigenvalue, 
Ao = 0, is given by 

|V#> = ^NW-*i|3 a >], Af =tl + tl (9) 

This is a CPT state that does not contain a contribu- 
tion from the intermediate state |2 Q ). The other two 



eigenstates 

m = -^=[h\i a ) - \±\2 a ) +t 2 \3 a )], 

Af± = t\ + )?± + t\ = 2Wo, 

with corresponding eigenvalues A± = ±\A? + contain 
all three states \j a }- If for a given coupling strengths t\ 
and t 2 the system is prepared in the CPT state ©, it 
will remain in this state as long as the couplings are con- 
stant in time. But even for time-dependent couplings, 
the system initially prepared in the CPT state can adi- 
abatically follow this state, provided the tunneling rates 
change slowly enough. More quantitatively, the nonadi- 
abatic coupling between the eigenstates of Hamiltonian 
([5]) is small, if during the evolution the transition am- 
plitude (V'±|i/'o) remains much smaller than the energy 
separation between the corresponding eigenstates Q, 

|(^|C)|«|A ± -A |. (10) 

Our objective is to transfer the electron from the first to 
the last QD using the time-dependent (pulsed) tunnel- 
couplings. From Eq. ^ one can see that if at an early 
time the tunnel coupling t 2 is switched on while t\ <C t 2 , 
the CPT state coincides with the initial state |1 Q ). One 
then slowly (adiabatically) decreases t 2 while increasing 
t\, so that at a later time t\ ^> t 2 and the CPT state 
coincides with the final state |3 a ). Assuming that t 2 
and ti are represented by partially overlapping pulses, 
each having temporal width t w , the adiabaticity condi- 
tion (UHl) requires ii n f x r w > 1. 

In the field of atomic/molecular physics, this tech- 
nique, involving the so-called counterintuitive sequence of 
pulses, is known as the stimulated Raman adiabatic pas- 
sage (STIRAP) that is commonly used for coherent pop- 
ulation transfer in three-state systems Q ■ We note that 
the solid-state implementations of the CPT and STIRAP 
in a pair of coupled quantum dots driven by two electro- 
magnetic fields has been proposed in [24[ . The single elec- 
tron transfer in a chain of three QDs via counterintuitive 
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FIG. 3: Time-evolution of a single-electron wavepacket in a chain of TV = 9 QDs with time-dependent, counterintuitive tunneling 
rates, (a) Population transfer is incomplete, |^1jv| 2 — 0.7, when the adiabatic condition is not very well satisfied, (b) Almost 
complete population transfer is achieved, | j4at | 2 — 0.97, when the adiabatic condition is better satisfied by doubling the temporal 
widths of the pulses and the total interaction time (note the different scales of the time axis in (a) and (b)). The insets show 
the time-dependence of even and odd tunneling rates and the populations of the first and last QDs. 



pulsing of tunnel-couplings as discussed above has been 
studied by Greentree et al. in where it was termed 
coherent tunneling by adiabatic passage (CTAP). These 
authors also considered the extension of CTAP to mul- 
tidot systems employing the so-called straddling scheme 
of Other schemes for adiabatic electron transport in 
tunnel-coupled QDs have been discussed in [25j]. 

Another extension of the STIRAP technique to sys- 
tems containing more than just three states has been 
given in [l(|. This scheme can easily be adapted to our 
system, as described below. We thus consider a chain of 
N sequentially coupled QDs and assume that the individ- 
ual tunnel couplings can selectively and independently be 
manipulated. When N is odd, i.e. N = 3, 5, 7, . . ., the 
Hamiltonian |5]| has a CPT eigenstate 



[£2*4 • ■ ■ ^N-l |1q) + (—1)^1^4 • ■ ■ tpf-l |3 Q ) 



+ ... + (-l) J t 1 t 3 ...t N - 2 \N a )}, 



(11) 



with eigenvalue Ao = 0. Thus the amplitude of the ini- 
tial state |1 Q ) is proportional to the product of all the 
even-numbered tunnel- couplings, while the amplitude of 
state \N a ) is given by the product of all odd- numbered 
tunnel-couplings, divided by the normalization parame- 
ter Af = (t 2 U . . . iAr_i) 2 + . . . + (iit 3 . . . ijv_ 2 ) 2 - There- 
fore, if all the even-numbered tunnel- couplings are pulsed 
together first, the CPT state (jTTJ) would coincide with the 
initial state |l a ). This is then followed by switching-on 
all the odd-numbered tunnel- couplings, while the even- 
numbered ones decrease, which will result in a complete 
transfer of electron wavepacket to the state \N a ). If we 
assume that these two families of pulses are described 
by common shape functions, t 2 , £4, . . . , i/v-i — ^oven and 



ii, £3, . . . , = ^odd, Eq. pip takes a compact form 

i i 

l^o) = -7Tr£H°<") n *^K 2n + 1 U ( 12 ) 

J 

KT ST* f 2n .2(J-n) 

• yv ° — ° dd cvcn ' 

71=0 

which makes the above discussion more transparent. In 
particular, complete population transfer from the initial 
state |1 Q ) to the final state \N a ) can be achieved by 
applying first the t cvcn pulses and then the i dd pulses, 
the two sets of pulses partially overlapping in time, as 
shown in Fig. [3] 

In order to minimize the nonadiabatic coupling of the 
CPT state to other eigenstates of the system, the rate of 
change of t QVCn and i ddj given approximately by the in- 
verse pulse- width t^ 1 , should be small compared to cor- 
responding eigenenergies |A| ~ |icvcn + iodd|, which yields 
the same condition as above, t™^£ odd T w ^> 1. One can 
see from the results in Fig.[3](a), which were obtained pre- 
cisely for this reason, that when this condition is not very 
well satisfied, the population transfer is incomplete. As 
expected, when the tunneling rates are pulsed for longer 
times, or, equivalently, have larger amplitudes, the adia- 
baticity condition is satisfied better, resulting in the com- 
plete population transfer from the initial to the final dot 
of the chain, as seen in Fig. [3£b). The remarkable advan- 
tage of this method over the one described in the previous 
Section is that as long as the two sets of partially overlap- 
ping pulses are strong enough, the adiabatic transfer of 
population is expected to be robust with respect to small 
uncertainties and fluctuations of tunneling rates, just like 
its atomic/molecular counterpart in Rcfs. [1, 0, G3 ■ On 
the other hand, the electron transfer via effective col- 
lective 7T pulse can be achieved with smaller tunneling 
rates and/or reduced interaction times, provided a pre- 
cise control of the tunneling amplitudes and timings is 
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possible. Depending on the characteristics of the partic- 
ular system, one or the other method may prove to be 
more practical. 

V. CONCLUSIONS 

In the above Sections, we have studied the dynamics 
of a single-electron transport in a linear array of tunnel 
coupled quantum dots. We have identified two regimes 
under which a complete coherent transfer of electron 
wavepacket between the two ends of the array can be 
achieved. Our results could be used for reliable infor- 
mation exchange between distant parts of an integrated 
quantum computer. As already noted in the Introduc- 
tion, one of the difficulties with the existing proposals 
for integrated QD based QCs [13, Ell is that the qubits 
(electron spins) interact with the nearest neighbors only, 
and there is no efficient way of transferring the infor- 
mation between distant qubits. As a way around such 
difficulties, one can envision an integrated quantum reg- 
ister composed of a large number of sub-registers, each 
containing two or more adjacent qubits, represented by 
spins of single electrons in individual QDs. The sub- 
registers are embedded in a two-dimensional array of 
empty QDs. As we have shown in an earlier publication 
[20] , through the mechanism of transient Heisenberg cou- 
pling, combined with the control of tunnel- coupling be- 
tween the dots studied in this paper, this two-dimensional 
grid could realize a flexible quantum channel, capable of 



connecting any pair of qubits within the register. Thus, 
to transfer the information, one connects distant sub- 
registers by a chain of empty QDs and applies one of the 
protocols described in the previous Sections to achieve a 
non-dispersive transfer of the qubit, followed by its con- 
trolled entanglement with a target qubit [13]. Note that 
this scheme is analogous to a proposal for an integrated 
ion trap based QC [26[, where, in order to circumvent the 
difficulties associated with a single large ion trap quan- 
tum register, it has been proposed to use many small 
sub-registers, each containing only a few ions, and con- 
nect these sub-registers to each other via controlled qubit 
(ion) transfer to the interaction region (entangler) repre- 
sented by yet another ion trap. 

We should note that the coherent electron dynamics in 
arrays of tunnel-coupled QDs bears many analogies with 
spin-wave dynamics in spin chains [ 27j or electromag- 
netic field dynamics in periodic photonic crystals = 2$, ;29t] , 
where some of the effects described above should be ob- 
servable. With an unprecedented control over system 
parameters, arrays of QDs doped with more than one 
electron allow for studies of numerous coherence and cor- 
relation effects in many-body physics. 
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